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We demonstrate how a topological atom laser can be realized by output coupling a trapped vortex state with a 
Raman scattering process. We find a linearized analytic solution from which a generalized resonance condition 
for Raman output coupling is developed. Using numerical simulations of a two component Gross-Pitaevskii 
equation in two and three dimensions the output beam from a trapped central vortex state is analyzed for cases 
of pulsed and continuous coupling where the vortex core is either transverse or parallel to the direction of 
propagation. We show how the parameters of the Raman light fields control the spatial phase of the output 
beam. 
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I. INTRODUCTION 

One of the most important technological developments to 
arise to date from dilute gas Bose-Einstein condensation is the 
atom laser Q10000- Analogously to its optical counter- 
part, the atom laser produces highly coherent, directed matter 
waves, and properties of this output beam, such as temporal 
and spatial coherence, and beam divergence have been char- 
acterized experimentally L£}|2l|g]- The typical atomic laser 
configuration uses a confining magnetic potential to act as a 
cavity for the laser mode, which is occupied by a Bose Ein- 
stein condensate. Radio-frequency or optical Raman transi- 
tions are then used to coherently transfer atoms into untrapped 
hyperfine states which can propagate freely away from the re- 
maining trapped atoms Jj§|. Usually the condensate is taken 
to be in the ground motional state of its confining potential, 
though it may possess thermal excitations fioll . 

In this paper we consider the matter wave output from a 
condensate in a topologically excited state. To be definite, we 
take the condensate to be in a central vortex state, though our 
formalism is sufficiently general to apply to a broader range 
of topological states including non-stationary states such as 
vortex arrays [e.g., see II 111 1. We present a T = OK Gross- 
Pitaevskii treatment of the topological atom laser dynamics, 
and give numerical results that characterize the behavior over 
a wide parameter regime, for two different scenarios: pulsed 
and continuous scattering of the matter field with the output 
beam direction transverse to the vortex core, and continuous 
scattering in the direction parallel to the vortex core. We ob- 
tain a linearized analytic solution for the output matter wave, 
and use it to determine the resonance condition for the Raman 
process and to characterize the phase properties of the atom 
laser. 



n. FORMALISM 

A. Raman Coupling 

In this paper we consider a Raman output coupling mech- 
anism of the type experimentally demonstrated by Hagley et 
al. 120. In that scheme the absorption and emission of opti- 



cal photons from far detuned laser fields coherently transfers 
atomic population between hyperfine states and imparts suffi- 
cient momentum for the output coupled atoms to recoil rela- 
tive to the trapped condensate. 

The ideal situation for an atom laser is to have only two 
hyperfine states coupled: (1) a weak field seeking state which 
we label |1) which is magnetically trapped by the confining 
potential into the laser modes; (2) an untrapped state |2) in 
which atoms can freely propagate to form a directed output 
matter wave. 

The states |1) and \2) have an energy difference E2 — E1 = 
UUJ21, and are coupled through an intermediate state |e) by two 
laser fields, with frequency u)\ and wavevector ki, and fre- 
quency L02 and wavevector k2 respectively. The atom makes a 
transition from state 1 1) to state |2) by absorbing an u)\ photon, 
and emitting an L02 photon, thus transferring kinetic energy 

Huj = %(u)\ — UJ2) ~ ftUJ21 (1) 

and momentum Kk = ft(ki — ^2) to the atom's center of 
mass motion. The laser fields are assumed to be sufficiently 
far detuned from the individual transitions that the effects of 
spontaneous emission can be ignored, and the laser detunings 
can be represented by a single value A. We can then adiabat- 
ically eliminate the intermediate state |e), and represent the 
laser field coupling between states |1) and |2) by a two pho- 
ton Rabi frequency V = ilifl^/^A, where fij is the single 
photon Rabi frequency associated with the ojj radiation field. 

The condensate atoms will be scattered into the output state 
only if the Raman energy transfer is approximately resonant 
i.e. uj w LUk where hujk — h 2 k 2 /2m is the recoil energy 
due to the momentum transfer to the atoms. We note that the 
quadratic dependence of the recoil energy on k makes subse- 
quent Raman transitions into other Zeeman sublevels nonres- 
onant. Thus provided V is small compared to w^, any sub- 
sequent coupling will be negligible, and the system can be 
regarded as being effectively two state. 

B. Matter Wave Evolution 

The condensate atoms in internal state |1) have center of 
mass wavefunction ipi(r) and are trapped in the cavity res- 
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FIG. 1: Raman couplings between internal Zeeman sublevels and 
momentum states in an F = 1 hyperfine manifold, appropriate for 
Rb 87 and Na 23 . The frequencies of the two Raman lasers are toi and 
u>2 respectively, and the other quantities are defined in the text. 



onator mode by a harmonic potential Vt (r). Correspondingly, 
the atoms in state |2) have wavefunction ^a(r). The stationary 
modes of the trap (or resonator) satisfy the Gross-Pitaevskii 
equation 



We have taken all scattering lengths between and within states 
to be degenerate, and we note that the normalization condition 
for the wavefunctions is J dr [\ipi\ 2 + IV^I 2 ] = 1. In this 
paper we solve Eqs. (0 and Q numerically in two and three 
spatial dimensions and also obtain an analytic solution for the 
output wave. 



III. LINEARIZED ANALYTIC SOLUTION 

We will consider the regime where the depletion of the 
trapped condensate is negligible on the time scale of the out- 
put coupling. This condition will allow us to consider a pertur- 
bative solution for the output beam amplitude, valid to first or- 
der in V. At this order, we can ignore depletion of the trapped 
condensate, and write ■0i(r,t) = yj p(r) exp (iSo(r) — ifit) 
where p(r)(=|-0i(r, 0)| 2 ) and Sq(t) are the density and phase 
profiles of the initial trapped state respectively. We will also 
assume that we can neglect the collisional interaction of atoms 
in state \2) on themselves, and on atoms in state |1) . This 
requires that the density of the scattered atoms is small com- 
pared to the trapped atom density in the region of overlap, 
which is consistent with the negligible depletion assumption. 
Transforming ip2 to an interaction picture (to remove the free 
particle-like motion) defined as ip2^,t) = ^2(1", t) exp(ik • 
r — iojt), Eq. @ can now be written in its linearized form as 
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= VVi +V T (r)^i +w|^i| 2 V'i, (2) 
im 

where Tip is the chemical potential and the collisional interac- 
tion strength w is given by 4irh 2 anTVo / m with an the s-wave 
scattering length for state |1) and No the number of conden- 
sate atoms. 

For a topological atom laser we require that the trapped 
condensate is in a topological state, whose properties will be 
transferred to the output matter wave. For simplicity we con- 
sider the case where the initial condition for the condensate 
is a solution of Eq. (|2jl that possesses a phase winding: in 
particular we will assume a central vortex state with a single 
quantum of circulation. 

The time evolution of the trapped condensate and the out- 
coupled matter fields at T = OK is given by the coupled 
Gross-Pitaevskii equations fl2ll 
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(4) 



where S = uj — cok is the Raman detuning from free particle 
resonance, K = Kk/m is the recoil velocity, and in accor- 
dance with the linearized treatment we have ignored nonlinear 
terms involving ?/>2 ■ 

If the velocity spread in the trapped condensate is small 
compared to K then we can ignore the effects of diffusion 
(about the central momentum frk) of the scattered wave packet 
by neglecting the Laplacian term in Eq. l|5}. Noting that the 
local velocity of the trapped condensate at position R is 



v(R) = 



W R So(R) 



(6) 



this approximation is valid where v(R) 2 <c K 2 . For typ- 
ical experimental parameters this condition will be well sat- 
isfied everywhere on the trapped condensate except close to 
the vortex core where the velocity field diverges. Making this 
approximation we obtain the formal solution for the output 
coupled wave packet 



where 



ds e i@( - r ' t ' s W{s)y/p(r + K(a - t)), (7) 



6(r,t,s) = So(r + K(s-t))-fj,s + (t-s)6 



~J ds' p(r + K(s' -<)). 



(8) 
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In previous work we have presented a similar approximate so- 
lution for Bragg scattering (where only a single internal state 
is involved) and have verified that it provides a good represen- 
tation of the behavior of the Gross-Pitaevskii equation over a 
wide parameter regime 1 13J] . 

The interpretation of Eq. is similar to the Bragg case. As 
the scattered packet moves across the trapped condensate, the 
amplitude ip2 at a given point (stationary in the frame mov- 
ing with velocity K ) is built up from the sum of contribu- 
tions coupled in from the trapped condensate at successive 
points along the trajectory R = r + K(s — t). The contri- 
bution coupled into ip2 at time s from position R, has moved 
to position r at time t and and has acquired a net phase of 
0(r, t, s). If the phase term 9 varies sufficiently slowly along 
a given trajectory, then the cumulative contributions of matter 
scattered from the trapped condensate will interfere construc- 
tively. Maximal scattering occurs when S is chosen so that 
there is a point of stationary phase along the trajectory (i.e. 
d<d(r, t, s)/ds = 0), which gives the generalized Raman res- 
onance condition 

5 « [V R S (R) • K - M + <MR)/% =r+K(s _ t) • (9) 

It should be stressed that Eq. (|9j is a local condition, which 
will apply only at certain positions in the condensate. This 
leads to the important property that the Raman scattering pro- 
cess can be spatially selective. The resonance condition Eq. 
differs from the Bragg resonance condition we previously 
derived 1131 . and gives rise to the following physical inter- 
pretation. Immediately after the Raman ejection process, an 
atom at R will have velocity K + v(R) [see Eq. (|6j]. The 
atom will also have center of mass energy %(u> + fi), arising 
from the chemical potential released from the trapped conden- 
sate (hp), together with an additional energy huo gained from 
the radiation field. This energy is divided between kinetic en- 
ergy m(K + v(R)) 2 /2 and potential energy wp(R) arising 
from interaction with the atoms remaining in the trapped con- 
densate. Energy conservation during the Raman process is 
therefore expressed as 

h(w + n) = y(K + v(R)) 2 +wp(R). (10) 

One can easily see that Eq. dlOt is equivalent to Eq. pro- 
vided we neglect the term in v(R) 2 , consistent with the valid- 
ity condition v(R) 2 <C K 2 used in deriving the analytic solu- 
tion. It is also convenient below to interpret the term VSo • K 
as the Doppler shift due to the local velocity of the condensate. 
Two particular cases of the resonance condition are worth ex- 
amining: 

(i) For a non-interacting ground state (i.e. So = 0, w — 
0) the resonance condition (|9} is 6 = — /i. This is in 
contrast to the result for Bragg scattering for this case, 
where resonance is at <5 = 0, with the difference due 
to the fact that the Raman process ejects the atom from 
the trapped initial state, so releasing the energy of that 
state. 

(ii) For an interacting ground state (i.e. So = 0) the res- 
onance condition is S = wp(R)/fr — p. By density 
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FIG. 2: Pulsed output matter wave from a 2D, m = — 1, central 
vortex state in a radially symmetric trap at t = 0.71to- (a) Temporal 
dependence of the Raman coupling V(t). (b) Output density profile, 
(c) Output phase profile. Raman output coupler parameters: k = 
30\iy/xo, uj = 925u;t, S = 25ujt, and Vm a x = 25ll>t- Condensate 
parameters: w — 500 wo and p, = 9.2ojt- Quantities are expressed 
in harmonic oscillator units with to = 1/wt, %o = yS/2mwr, 
and wo = TlujtXq, where lot is the harmonic trap frequency for the 
trapped condensate. 

weighting the contributions from all R across the con- 
densate (assuming a Thomas-Fermi profile) we find that 
the maximum amount of atoms will be scattered when 
<5 = -3/x/7. 



IV. PULSED OUTPUT COUPLING 

Here we consider a configuration of Raman fields which 
are pulsed in time with a scattering direction (i.e. K) per- 
pendicular to the vortex core. With this choice, one axis of 
the plane containing the characteristic spatial and phase vari- 
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ations of the vortex is parallel to the direction of the beam 
propagation. For convenience we present 2D simulations of 
this geometry with the Raman momentum difference vector 
along the y-direction, i.e. k = fey. The initial trapped state is 
assumed to be of the form i/ji(x, y) — R(r)e 1 ^, where r and 
4> are the radial coordinate and polar angle respectively. 

In Fig. [2]we show the output matter wave (i.e. ^)%) after the 
application of four Raman pulses. The length of each pulse t p 
is chosen to be sufficiently short that the atoms hardly move 
during the pulse, i.e. r p <c R c / K (where R c is the condensate 
radius), while the time between pulses is sufficiently long that 
the scattered atoms can traverse the trapped condensate before 
the next pulse is applied. The Raman intensity envelope V(t) 
is shown in Fig. EJa), and the resulting output coupled den- 
sity and phase are shown in Figs. |2jb) and (c). It is clear that 
this procedure produces copies of the trapped vortex state in 
the output state, including a 2ir phase circulation, a result that 
we confirm analytically below. We note that because {p2 is 
defined in an interaction picture with respect to the free parti- 
cle center-of-mass motion, the rapid phase gradient due to the 
mean wave packet velocity K has been removed, allowing the 
phase information in Fig. |2jc) to be more easily discerned. 
For this reason all the phase plots we show in this paper will 
be given for the interaction wavefunction. 

The Fourier frequency width of each short pulse in our sim- 
ulation is ~ 210 Ut, which is much larger than the typical fre- 
quencies associated with meanfield and Doppler shifts. In this 
transient temporal regime, the resonance condition Eq. i|9} is 
inapplicable, however the frequency spread is still sufficiently 
narrow that there is no significant coupling into other mag- 
netic sublevels. As long as the central frequency is close to 
resonance (u w u>k) this temporal-limited frequency spread 
causes the Raman coupling to be resonant everywhere on the 
trapped condensate facilitating the complete copying of the 
wavefunction into the output component. The analytic solu- 
tion Eq. contains the vortex output result of Fig. [2] as we 
can show by modeling an individual Raman pulse at time tj 
as a delta function Vo6(t — tj). For simplicity we shall choose 
Vq to be real, and then using Eq. Q we obtain for t > tj (in 
the non-interaction picture) 

V> 2 (r,i) = ~V oy /p{r - TLj) exp(ik • r) exp(»0) (11) 

where R 3 = K(t — tj), and 

9 = -ut + S (r - Rj) - fdj + (t - tj)S 

™ / ds' p(r + K(s' -t)). (12) 
ft Jtj 

The density distribution of the output field in Eq. il li is the 
same as the initial vortex state, but displaced by Rj (and mul- 
tiplied by Vq/2 ). The phase 9 of the field can be written 

= [So(r-Rj) -fdj] -ojtj -u k (t-tj) 

-- / ds'p(r + K(s' - 1)) (13) 
ft J t] 

which has the following interpretation. The term in square 
brackets represents the phase of the trapped vortex at time tj , 
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FIG. 3: Continuously output coupled matter wave from a 2D, m = 
— 1, central vortex state of a radially symmetric trap at t — 0.71to- 
with continuous Raman coupling, (a) Density profile. Dashed line 
indicates steady state beam region (see text). The solid lines indicate 
density contours of the trapped vortex state, (b) Phase profile. Raman 
parameters: V — 2, ui — 925lut, k = 30k H /:co and S = 25u>t- 
Condensate parameters: w = 500it)o and fi = 9.2wt- 



which is transferred to the output vortex along with the phase 
of the Raman coupling (u)tj) at the time of the pulse. In the 
absence of collisional interactions the vortex would behave 
as a free particle of momentum k, and energy LOk (we have 
neglected diffusive kinetic energy) and its phase would incre- 
ment as LOk (t — tj) . The integral term in Eq. Jl 3i represents 
the correction to the free particle evolution during the time the 
output field is transiting the trapped vortex. Notice that the 
integrand cuts off once |r + K(,s' — t) \ > R c ,so that this con- 
tribution to the phase becomes constant once the output vortex 
has separated from the initial vortex. The rotation of the line 
of zero phase between successive vortices in Fig. |2jc) is thus 
described by the first 3 terms of Eq. Jl 31 . 
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V. CONTINUOUS OUTPUT COUPLING 



In Fig. [3] we show the output matter wave following the 
application of a long duration, low intensity Raman pulse to a 
vortex state - which we will refer to here as continuous output 
coupling. This result displays typical behavior of the scattered 
solution in the continuous regime: the matter wave is prefer- 
entially scattered from a selected spatial region of the trapped 
condensate and has an asymmetric density distribution trans- 
verse to the direction along which it propagates. The shape 
of the output density profile depends on the detuning S, and 
by suitably adjusting this parameter, matter can be selectively 
coupled out from either side of the trapped vortex. Here the 
frequency width of the finite duration Raman pulse is suffi- 
ciently narrow that the generalized resonance condition l|9) is 
applicable. From this resonance condition, the spatial asym- 
metry of the output density profile can be understood as aris- 
ing from the Doppler term VrS'o(R) • K. A unit circulation 
central vortex state with angular momentum number m = ±1 
has a phase profile of the form So = ±0, where tf> is the az- 
imuthal angle. The Doppler effect gives a local shift of the 
resonant frequency; on the side of the vortex where the cur- 
rent flow runs parallel to the scattering direction the resonant 
frequency is shifted upward, while on the other side the flow 
is anti-parallel and the frequency is shifted downward. By 
choosing S above or below the central resonance value, it is 
possible to selectively scatter from either side of the vortex. 

To investigate the properties of the continuously output 
matter wave in more detail we consider the steady state re- 
gion of the output coupled matter wave. This region, which 
we indicate pictorially in Fig. 13a), is where the output beam 
has constant transverse density profile, that we shall refer to 
as the steady state beam profile. From the formal solution of 
Eq. Q it can be shown that for two points r and r + otK, 
both within the steady state beam region, then the values of 
the wavefunction at these points are related by 



^(r + aK,f)=^(r,f)e l(i+ " )a . 



(14) 



This result indicates that the steady-state velocity of the output 
matter wave is K + Vq , where we define 



V = A^k. 

771 K 



(15) 



Remarkably, this is uniform across the transverse profile, a 
result we can understand using Eq. which shows that at 
the position where the Raman transfer takes place, the trapped 
condensate has a local velocity component parallel to K given 
by 
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FIG. 4: Population scattered into an output matter wave from a 
vortex after excitation by a continuous Raman pulse of duration 
t — 0.71io- Raman parameters: k = 30\i y /xo and V = 2ut- 
Condensate parameters: w = 500^0 and fi = 9.2wr- 



moved away from the condensate, this potential energy is con- 
verted to kinetic energy, consistent with an atom of velocity 
K + "Vq. We note that the transverse components of the lo- 
cal ejection velocity cancel out in the final scattered wave, 
due to the symmetry of the vortex velocity distribution. For 
a general straight-line trajectory across the vortex, there are 
two positions of local resonance with the same density p(R), 
but with opposite transverse local velocities. It is also worth 
emphasizing that the control available over the output mat- 
ter wave velocity is in contrast with the case of optical lasers, 
where the group velocity of the light cannot be altered without 
changing the properties of medium in which it propagates. 

We have numerically verified that Eq. fl!4i is in good 
agreement with numerical simulations of the Gross-Pitaevskii 
Equation over a wide regime. We note that changing the Ra- 
man detuning S to control the phase gradient of the output 
coupled beam also affects the efficiency at which the atoms 
are scattered, and to characterize this we show the scattered 
population, i.e. n,2 = J G?r|-02(r, t)\ 2 in Fig. |4] This result 
shows that in this case the frequency width of the Raman cou- 
pling (from S ~ — SQlut to 30cjt) is sufficient to allow the 
output velocity to be changed over a range of 1.0xqu>t with- 
out significant attenuation of the output matter wave. In these 
results the Doppler effect is the dominant broadening mecha- 
nism, so that the approximate width of the Raman transition, 
Au, will be 



Au = 



kAp 



(17) 



v(R) 



h (6 



fi - wp(R)) £ 



K 



(16) 



The velocity of the atom (in direction K) immediately after 
the Raman transfer is therefore K + v(R)y. We showed in 
section III that this local ejection velocity gives a kinetic en- 
ergy consistent with an atom residing in a potential of ijjp(R) , 
(due to the trapped condensate atoms). Once the atom has 



where Ap is the momentum width of the condensate in the 
scattering direction. For the results in Fig. |4]Eq. illl gives 
the full width at half maximum Aoj = 60.8ujt, where Ap is 
numerically evaluated from the initial vortex eigenstate (the 
initial state for ipi). 

Additionally, it should be noted that the scattered popula- 
tion (712) is a convenient experimental observable for mea- 
suring the linear response of the condensate to the perturbing 
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FIG. 5: Output matter wave at t = 1.4ms from a 2.2 x 10 5 atom 
23 Na condensate in a m — 1 central vortex state with /i — 1.15kHz 
in a spherically symmetric trap of frequency 80Hz. Raman output 
coupler corresponds to 589nm laser fields at 120° with S = and 
a two-photon Raman frequency of V — 2ir x 480s -1 . (a) 3.4 x 
10 12 atom/cm 3 density iso-surface of the output matter wave, (b) 
and n phase iso-surfaces of the output matter wave. 

Raman potential. This technique for probing the condensate 
is similar to Bragg spectroscopy, which was first used on con- 
densates by Stenger et al. | 14], except that for the Raman case 
the condensate is scattered into a different internal state. This 
technique has several advantages over Bragg spectroscopy for 
conducting high precision measurements of condensate prop- 
erties which we will explore elsewhere. 

Finally we consider continuous output coupling of a central 
vortex in a direction parallel to the vortex core. Because this 
arrangement has the superfiuid current flow perpendicular to 
K, the Doppler term and hence the output coupling is not spa- 
tially selective 1 17]. This means that the vortex profile will be 
reproduced as the transverse profile of the atom laser output, 
so that the atom laser will have a similar intensity and phase 
profile to that of a TEMqj mode of an optical laser, commonly 



known as the donut mode lll6ll . 

In Figs. |5] (a) and (b) we show phase and density iso- 
surfaces of a vortex output-coupled in the manner described 
above, for realistic experimental parameters. The noticeable 
beam divergence in the stable beam region, which extends 
from about z = lOfim to z — 60^m, is enhanced by the cen- 
trifugal forces of the rotating condensate beyond that expected 
simply from repulsive wave packet spreading 

The helicity of the phase iso-surface arises from the addi- 
tion of the phase gradient along the propagation direction (due 
to the Raman output coupling) to the vortex-like phase profile 
the matter wave has in the transverse direction. The phase 
gradient of the matter wave and hence the phase helicity is 
controlled by Raman detuning, as given in Eq. (I14> . We have 
numerically verified that for a case with the same parameters 
as in Fig. [5] except using 6 ~ — /i, the helicity disappears 
as the phase iso-surfaces are approximately vertical. Also for 
5 < — /j, the helicity changes sign corresponding to the phase 
iso-surfaces spiraling in the opposite sense. 

The helical structure of the phase could be investigated by 
an interference experiment, such as superimposing the topo- 
logical atom laser beam with a co-propagating plane phase 
atom laser. 



VI. CONCLUSION 

In this paper we have introduced the idea of a topological 
atom laser and considered its behavior using numerical and 
analytic approaches over a broad parameter regime. We have 
developed a linearized solution for the output coupled matter 
wave from which we have determined a spatially dependent 
resonance condition for Raman scattering and characterized 
how the output beam phase properties relate to the Raman de- 
tuning and the condensate chemical potential. 



VII. ACKNOWLEDGMENTS 

This work was supported by the US Office of Naval Re- 
search, the Advanced Research and Development Activity. 
PBB and RJB acknowledge support from the Marsden Fund 
of New Zealand under contracts PVT902 and PVT202. 



[1] M.-O. Mewes, M. Andrews, D. Kurn, D. Durfee, C. Townsend, 
and W. Ketterle, Phys. Rev. Lett. 78, 582 (1997). 

[2] B. Anderson and M. Kasevich, Science 282, 1686 (1998). 

[3] E. Hagley, L. Deng, M. Kozuma, J. Wen, K. Helmerson, S. Rol- 
ston, and W. Phillips, Science 283, 1706 (1999). 

[4] I. Bloch, T. W. Hiinsch, and T. Esslinger, Phys. Rev. Lett. 82, 
3008 (1999). 

[5] J. Martin, C. McKenzie, N. Thomas, J. Sharpe, D. Warring- 
ton, P. Manson, W. Sandle, and A. Wilson, J. Phys. B 32, 3065 
(1999). 

[6] M. Kohl, T. Hansch, and T. Esslinger, Phys. Rev. Lett. 87, 



160404 (2001). 

[7] M. Trippenbach, Y. Band, M. Edwards, M. Doery, P. Julienne, 

E. Hagley, L. Deng, M. Kozuma, K. Helmerson, S. Rolston, 

et al., J. Phys. B 33, 47 (2000). 
[8] Y. L. Coq, J. H. Thywissen, S. A. Rangwala, F. Gerbier, 

S. Richard, G. Delannoy, P. Bouyer, and A. Aspect, Phys. Rev. 

Lett. 87, 170403 (2001). 
[9] R. Ballagh and CM. Savage, in Proceedings of the Thirteenth 

Physics Summer School: Bose-Einstein Condensation: Atomic 

Physics to Quantum Fluids, edited by C. Savage and M. Das 

(World Scientic, Singapore, 2000). 



7 



[10] S. Choi, Y. Japha, and K. Burnett, Phys. Rev. A 61, 063606 
(2000). 

[11] K. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. 

Rev. Lett. 84, 806 (2000). 
[12] M. Edwards, D. A. Griggs, P. L. Holman, C. W. Clark, S. Rol- 

ston, and W. Phillips, J. Phys. B 32, 2935 (1999). 
[13] P. B. Blakie and R. J. Ballagh, Phys. Rev. Lett. 86, 3930 (2001). 
[14] J. Stenger, S. Inouye, A. Chikkatur, D. Stamper- Kurn, 

D. Pritchard, and W. Ketterle, Phys. Rev. Lett. 82, 4569 (1999). 



[15] P. Blakie and C. W. Clark (2002). 

[16] A. Siegman, Lasers (University Science Books, 1986). 

[17] There is a residual resonance shift due to the inhomogeneous 
meanfield, but we have found that typically this does not give 
rise to significant spatial selectivity if the condensate chemical 
potential is of order or smaller in size that the Doppler width 
(see Eq. O)- 



